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We study dynamical supersymmetry breaking and the transition point by non-perturbative lattice 
techniques in a class of two-dimensional N = I Wess-Zumino model. The method is based on 
the calculation of rigorous lower bounds on the ground state energy density in the infinite-lattice 
limit. Such bounds are useful in the discussion of supersymmetry phase transition. The transition 
point is determined with this method and then compared with recent results based on large-scale 
Green Function Monte Carlo simulations with good agreement. 
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1. Introduction 

An important issue in the study of supersymmetric models is the occurrence of non-perturbative 
dynamical supersymmetry breaking l . The problem can be studied in the N = I Wess-Zumino 
model that does not involve gauge fields and is thus a simple theoretical laboratory. Since the 
breaking of supersymmetry is closely related to the symmetry properties of the ground state, we 
will adopt a Hamiltonian formulation of the model. 

Let us remind the (continuum) N = \ super algebra, {Q a ,Qp} = 2(/ 5 C) aj g. Since Pj are not 
conserved on the lattice, the super algebra is explicitly broken by the lattice discretization. A 
very important advantage of the Hamiltonian formulation is the possibility of conserving exactly 
a subset of the full super algebra [pp. Specializing to 1 + 1 dimensions, in a Majorana basis y = 
C = a 2 , 7i = ia-s, the algebra becomes: Q\ = Qj = P° = H and {Qi,Qi} = 2P l = IP. On the 
lattice, since H is conserved but P is not, we can pick up one of the supercharges, say, Q\, build 
a discretized version <2l and define the lattice Hamiltonian to be H = Q\. Notice that Q\ = H 
is enough to guarantee Eq > 0. The explicit lattice model is built by considering a spatial lattice 
with L sites. On each site we place a real scalar field cp n together with its conjugate momentum p n 
such that [p n ,(Pm] = —i$nsn- ^ ne associated fermion is a Majorana fermion \\f a n with a = 1,2 and 

{Ya,n, Yb,m} = 8a,b8n,m, Wa,n = Ya,n- 

The continuum 2-dimensional Wess-Zumino model is defined by the supersymmetric genera- 
tors involving the supeipotential V((p), 



Ql,2 = J dx 



(1.1) 



where cp(x) is a real scalar field and y(x) is a Majorana fermion. The discretized supercharge 
is IJ 

<2L=£ PnYhn ~ ( — ~ <P "' 1 +V(<Pn)) Yl,n (1-2) 

and the Hamiltonian takes the form 

2 



^ n=\ 



n 2 + ^+1^-1 + _ {xhn+i + h c ) + (_l)V(^) [2xlXn ~ 1) 



(1.3) 

where we replace the two Majorana fermion operators with a single Dirac operator % satisfying 
canonical anticommutation rules. 

The problem of predicting the pattern of supersymmetry breaking is not easy. In principle, 
the form of V(q>(x)) is enough to determine whether supersymmetry is broken or not. At least at 
tree level supersymmetry is broken if and only if V has no zeros. The Witten index [Q] can help 
in the analysis: If V (<p) has an odd number of zeroes then 1^0 and supersymmetry is unbroken. 
If V(q>) has an even number of zeroes, when / = we can not conclude anything. An alternative 
non-perturbative analysis for the case / = is thus welcome. The simplest way to analyze the 
pattern of supersymmetry breaking for a given V is to compute the ground state energy Eq through 



1 See [jj|] for recent reviews and a complete list of references. 
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numerical simulations and/or strong coupling expansion. On the lattice, accurate numerical results 
are available []|, §], although a clean determination of the supersymmetry breaking transition re- 
mains rather elusive. All the predictions for the model with cubic prepotential, V = (p 3 , indicated 
unbroken supersymmetry. Dynamical supersymmetry breaking in the model with quadratic prepo- 
tential V = fatp 2 + Ao was studied performing numerical simulations []|] along a line of constant A2, 
confirming the existence of two phases: a phase of broken supersymmetry with unbroken discrete 
Z2 at high Ao and a phase of unbroken supersymmetry with broken Z2 at low Ao, separated by a 
single phase transition. 

On the other hand, from the strong coupling analysis what comes out is the following: for odd 
q, strong coupling and weak coupling expansion results agree and supersymmetry is expected to 
be unbroken [Q]. This conclusion gains further support from the non vanishing value of the Witten 
index [Q]. For even q in strong coupling, the ground state has a positive energy density also for 
L — > 00 and supersymmetry appears to be broken. In particular, for V = A29 2 + Ao, weak coupling 
predicts unbroken supersymmetry for Ao < 0, whereas strong coupling prediction gives broken 
supersymmetry for all Ao- 

2. Numerical Simulations and Discussion 

We used two different approaches to investigate the pattern of dynamical supersymmetry 
breaking. In the first one, [0, H], the numerical simulations were performed using the Green Func- 
tion Monte Carlo (GFMC) algorithm and strong coupling expansion. The GFMC is a method 
that computes a numerical representation of the ground state wave function on a finite lattice with 
L sites in terms of the states carried by an ensemble of K walkers. Numerical results using the 
GFMC algorithm for the odd prepotential confirm unbroken supersymmetry. 

A more interesting case is the even prepotential. When V = Xiqr + Ao and for fixed A2 = 0.5, 
we may expect (in the L — > °o limit) a phase transition at Ao = kff^ (A2) separating a phase of broken 
supersymmetry and unbroken Z2 (high Ao) from a phase of unbroken supersymmetry and broken 
Z2 (low Ao). 

The usual technique for the study of a phase transition is the crossing method applied to the 
Binder cumulant, B. The crossing method consists in plotting B vs. Ao for several values of L. The 
crossing point Xq(L\,Lq), determined by the condition B{Xq,L\) = B(Xq,L2) is an estimator of 
Aq . The value obtained is showed in Fig. [j] and corresponds to Aq = —0.48 ±0.01 The main 
source of systematic errors in this method is the need to extrapolate to infinity both K and L. For 
this reason, an independent method to test the numerical results of ^ is welcome. 

The second method is based on the calculation of rigorous lower bounds on the ground state 
energy density in the infinite-lattice limit [^, Such bounds are useful in the discussion of super- 
symmetry breaking as follows: The lattice version of the Wess-Zumino model conserves enough 
supersymmetry to prove that the ground state has a non negative energy density p > 0, as its con- 
tinuum limit. Moreover the ground state is supersymmetric if and only if p = 0, whereas it breaks 
(dynamically) supersymmetry if p > 0. Therefore, if an exact positive lower bound Plb is found 
with < Plb < p , we can claim that supersymmetry is broken. 

The idea is to construct a sequence p ^ of exact lower bounds representing the ground state 
energy densities of modified lattice Hamiltonians describing a cluster of L sites and converging to 
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Figure 1: The Binder cumulant B vs. Ao for K = 200 and K = 500. Here Xi c) = -0.48 ± 0.01 . 



p{L) —> p i n the limit L — > oo. The bounds p( L ) can be computed numerically on a finite lattice with 
L sites. The relevant quantity for our analysis is the ground state energy density p evaluated on the 
infinite lattice limit p = lim^^ MS. it can be used to tell between the two phases of the model: 
supersymmetric with p = or broken with p > 0. 

In Ref . [^] we presented how to build a sequence of bounds p which are the ground state 
energy density of the Hamiltonian H with modified couplings on a cluster of L sites: given a 
translation-invariant Hamiltonian H on a regular lattice it is possible to obtain a lower bounds 
on its ground state energy density from a cluster decomposition of H, i.e., given a suitable finite 
sublattice A, it is possible to introduce a modified Hamiltonian H restricted to A such that its energy 
density Pa bounds p from below. The difference between H and H amounts to a simple rescaling 
of its coupling constants. The only restriction on H being that its interactions must have a finite 
range [§]. 

We compute numerically p at various values of the cluster L: if we find p( L ' > Ofor some 
L we conclude that we are in the broken phase. We know that p( L ' — s> p for L — > oo and the study of 
p( L ) as a function of both L and the coupling constants permit the identification of the phase in all 
cases. The calculation of p( L ) is numerically feasible because it requires to determine the ground 
state energy of a Hamiltonian quite similar to H and defined on a finite lattice with L sites. 

To test the effectiveness of the proposed bound and its relevance to the problem of locating 
the supersymmetry transition in the Wess-Zumino model we study in detail the case of a quadratic 
prepotential V = facp 2 + Xq at a fixed value A2 = 0.5 [^J. An argument by Witten |4j] suggest the 
existence of a negative number Aq such that p(Ao) is positive when Ao > A ( , and it vanishes for 
Ao < Aq*. Aq* is the value of Ao in which dynamical supersymmetry breaking occurs. In Fig. || 
we show a qualitative pattern of the curves representing p( L '(Ao). We see that a single zero is 
expected in p^(Ao) at some Ao = Ao(L). Since lim^cop^ = p, we expect that Ao(L) — * Aq for 
L — > 00 allowing for a determination of the critical coupling Aq . The continuum limit of the model 
is obtained by following a Renormalization Group trajectory that, in particular, requires the limit 
h -+ [§. 

The properties of the bound p( L )(Ao) guarantee that for L large enough it must have a single 
zero Aq (L) converging to Ao as L — » 00. In any case for each L we can claim that Aq > Aq (L). To 
obtain the numerical estimate of p( L )(Ao) we used the world line path integral (WLPI) algorithm. 
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Figure 2: Qualitative plot of the functions p(Ao) and p( L )(Ao). 



Cluster size L - 14 Cluster size L - 18 




Figure 3 : Plot of the energy lower bound p W ( j3 , T ) at L = 14 and L = 1 8 . 

The WLPI algorithm computes numerically the quantity p( L \p,T) = l ^ 1 "^ ( e T ' e T ~) I where 

the Hamiltonian for a cluster of L sites is written as H = Hi +H2, by separating in a convenient way 
the various bosonic and fermionic operators in the subhamiltonians H\ and Hi 2 - The desired lower 
bound is obtained by the double extrapolation p^ = lim«^ M limr-woP^(j3,r), with polynomial 
convergence ~ l/T in T and exponential in j8. Numerically, we determined p^ L \fi,T) for various 
values of j8 and T and a set of Ao that should include the transition point, at least according to the 
GFMC results. In Fig. [3] we plot the function p( L \fi,T) for the cluster sizes L = 14 and L = 18, 
various j8 and T = 50, 100, 150. Here we see that the energy lower bound behaves as expected: it 
is positive around Ao = and decreases as Ao moves to the left. At a certain unique point Aq (L), the 
bound vanishes and remains negative for Aq < Aq (L). This means that supersymmetry breaking can 

2 we do not report Hi and H2 here, see Ref. WU for details. 
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Figure 4: Plot of Xq(L) vs. 1 /L for L = 6, 10, 14, 18. The best fit with a quadratic polynomial in 1 /L gives 
Xq = -0.49 ± 0.06 that should compared with the best GFMC result obtained with K = 500 walkers. 

be excluded for Ao > min L A ( ) (L). Also, consistency of the bound means that Aq (L) must converge 
to the infinite-volume critical point as L — > oo. Since the difference between the exact Hamiltonian 
and the one used to derive the bound is &(l/L), we can fit Aq (L) with a polynomial in 1/L. This 
is shown in Fig. ^| where we also show the GFMC result. The best fit with a parabolic function 
gives Aq = —0.49 ±0.06 |^] quite in agreement with the previous Aq GFMC = —0.48 ±0.01 & In 
conclusion, both methods reported here are quite in agreement and confirm the existence of two 
phases separated by a single phase transition at Ao for the quadratic prepotential. 

Acknowledgments: The calculations have been done on a PC Cluster at the Department of 
Physics of the University of Parma. 
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